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A technique developed to measure the dissolution rate of glass fibers in the laboratory under conditions approximating those thought to 
exist in the lung has been applied to 30 glass compositions in the SiOz-AlzOj-BaOj-CaO-MgO-BaO-NazO system. Dissolution rates 
show an Arrhcnius dependence on temperature with an activation energy near 64 kJ/moL The rate increases with the pH value. Glass fiber 
dissolution is a leaching process which leaves behind a mechanically-weak layer consisting of Si02 and AI2O3. Dissolution morphology 
changes with dissolution rate, suggesting a change in dissolution mechanism. Dissolution rates vary with composition from 0.9 to 
887 ng^(cm^ h). AljOj strongly decreases the dissolution rate, BjOs, BaO, CaO, MgO, and NaaO increase the rate to about the same degree, 
and Si02 has little effect. 

Auflosung von Glasfasem in einer physiologischen Salzidsung 

In Anlehnung an Bedingungen, deren Vorliegen in der Lunge man annimmt, wurde ein Vcrfahren zur Auflosung von Glasfasem entwickelt 
und auf 30 Glaszusammensetzungen des Systems SiOz-AliOj-BaOa-CaO-MgO-BaO-NaiO angewendet. Die Aunosungsgeschwin- 
digkeiten zeigen einc Atrhenius-Abhangigkeit von der Temperatur mil einer Akti vicrungscncrgie von etwa 64 kJ/moL Die Geschwindigkeit 
nimmt mit dem pH-Wert zu. Die Glasfascrauflosung ist cin AuslaugprozeB, der eine mechanisch instabile Schicht aus Si02 und AI2O3 
zurucklaBt. Die bei der Auflosung an der Glasobcrflache sich ausbildende Morphologic hangt von der Auflosungsgeschwindigkeit ab, was 
die Annahme nahelegt, daB sich der Auflosungsmechanismus andert. Die Auflosungsgeschwindigkeiten variieren mit der Zusammen- 
setzung von 0,9 bis 887 ng/(cm^ h). AI2O3 reduziert die Auflosungsgeschwindigkeit dcutlich, wahrend B^Oj, BaO, CaO, MgO und Na/) 
diese in etwa gleichem MaBe beschleunigen, SiOi zeigt nur cinen gcringcn EinfluB. 



1. Introduction 

In recent years there has been increasing interest in 
the health effects of inhaled fibrous material. Asbes- 
tos fibers have been shown to cause cancer and other 
lung diseases, and there is concern that other fibrous 
materials may cause similar health problems as a 
result of their fibrous geometry. It has been hy- 
pothesized that the residence time of fibers in the 
lung has a strong influence on their ability to cause 
disease, and that this residence time depends strongly 
on the clearance rate and the dissolution rate of the 
fibrous material. For this reason, a number of studies 
have been conducted to measure the dissolution rates 
of various inorganic fibers under conditions thought 
to be similar to those in the lung [1 to 6]. These 
studies have surveyed many types of inorganic fibers. 
They suggest that the chemical composition of fibers 
can have a strong influence on the dissolution rate, 
but variations in the physical characteristics of the 
sample limit their usefulness in understanding the 
relation between chemical composition and dissolu- 
tion rate. Nor has there been any detailed study of the 
dissolution mechanism. This report presents the 
results of a systematic study of glass fiber dissolution 
rate and mechanism for glasses in the 
SiOi-AlsOa-BiOa-CaO-MgO-BaO-NaiO sys- 
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tern using glass fiber samples melted and formed in 
the laboratory to approximately the same uniform 
diameter. 

2. Experimental procedure and rationale 

Glasses were prepared from reagent grade materials 
by melting in a covered Pt-Rh alloy crucible for 3 h at 
1450 in an electric furnace. After each hour, the 
melt was removed from the furnace and stirred briefly 
to improve glass homogeneity. Fibers of uniform 
diameter were formed by standard techniques from a 
one-hole Pt-Rh alloy bushing. The fibers were cut 
into approximately 1 cm lengths and stored in sealed 
plastic containers in a desiccator. 

Fiber diameter measurements were made on 
polished cross sections of expoxy-embedded fibers 
using a Scanning Electron Microscope (SEM). The 
number of measurements made for each sample was 
generally in the range 20 to 50. Fibers near 10 txm in 
diameter were routinely used in the experiments due 
to ease of formation and handling, but for four 
compositions; fibers near 2 |im in diameter were also 
formed to apply the dissolution technique to fibers of 
respirable diameter and to confirm that the dissolu- 
tion-rate constant measured by this technique is a 
fundamental glass parameter independent of fiber 
geometry. 
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Rber dissolution has generally been measured by 
monitoring either fiber weight loss or the dissolution 
fluid composition. Neither method is ideal. The 
dissolution process at physiological pH values in- 
volves both selective leaching and network dissolu- 
tion. Fiber weight loss measurements do not allow 
distinguishing these two processes; complete analysis 
of the fluid for dissolved components is difflcult or 
impossible because some are already present in the 
fluid in concentrations far exceeding the amounts 
dissolved. Neither method addresses morphological 
changes, which may lead to physical breakup of the 
fibers. As Scholze [5] points out in a general review of 
the fiber dissolution literature, it is clearly desirable 
to monitor all possible changes: solution composi- 
tion, weight loss, chemical composition, and mor- 
phology of the fiber. For the present study, which 
aims at a general understanding of the effects of glass 
composition on dissolution, fiber weight loss data on 
many fiber compositions are sufficient to determine 
the basic compositional dependence. For some 
samples solution composition and dianges in fiber 
morphology were also measured in addition to weight 
loss. 

To prepare a sample for dissolution rate meas- 
urement, approximately 1.5 g of fiber was dispersed 
in deionized water by a 15 s treatment in a blender. 
The suspension was filtered through a 0.4 (un 
polycarbonate membrane filter backed by a plastic 
filter support pad in the bottom of a 37 mm diameter, 
three-section air-monitoring cassette. The cassette 
was reassembled, enclosing the sample between two 
0.4 (im polycarbonate filters. The mat of fibers 
generally filled the 6.6 cm^ space between the two 
filters. This procedure eliminates fiber bundles that 
could trap or exclude solution and locally alter the 
dissolution process. It allows the initial sample weight 
to be determined by difference fi*om the preweighed 
cassette and filters. And it keeps the fibers enclosed 
throughout subsequent procedures. After filtration, 
the cassette was dried for 2 d under vacuum at 50 
and then to constant weight in a glove box with 
CaS04 desiccant at room temperature. After final 
weighing, the sample was removed from the cassette 
and stored in a sealed plastic bag until needed for 
morphological studies. 

The composition of the fluid-contacting fibers in 
the lung is not known although some constraints can 
be put on it [7]. Most of the simulated lung fluid 
compositions that have been used in the literature are 
based on the composition of interstitial fluid as 
reported by Gamble [8]. That used in this study is a 
modification of Gamblers solution developed by 
Kanapilly [9], in which the calcium concentration has 
been reduced to avoid precipitation of calcium salts 
(table 1). To eliminate bacterial and algal growth, 
0.1 % NaNs has been added to this solution. Several 
experiments were run in different fluids to give some 
indication of the effect of solution composition. 



Table 1. Ptiysiologica] saline solution compositions in mg/I of 
solution and dissolution-rate constants for selected glasses 





moQifieo 


niuuiiicu 


moQitieu 






from F41 


piasiua ) 




7.9 


69.6 


95 








166 


Mg2+ 

Na"*" 




21.3 


25 


3363.3 


3516.8 


4267 




180.4 




\ 


a- 


4481.4 


4085.0 


3651 


HCOJ 


1647.3 


1963.3 


1216 


poj- 


114.2 


87.9 


287 


soj- 


48.0 


53.4 


32 




64.6 


64.6 




citrate 


37.9 


119.6 


3171 


glydne 


450.0 


118.0 




tartrate 




115.8 




lactate 




139.1 


107 


pyruvate 




136.1 




dextrin 






3600 


total protein . 


(4S0.0) 


(118.0) 


(> 72000) 


k in ng/(cm^h) 








glass no. 5 


27.4 


51.6 


25.4 


glass no. 22 


254.6 




144.1 



^) Outdated plasma from the American Red Cross, which has the 
following additives to whole blood: 3.7 g/t sodium citrate, 3.6 g/1 
dextrin, 460 mg/l citric acid, and 310 mg/1 NaH2P04. This plasma 
composition is based on typical values of human whole blood, 
plasma, or serum [10] with these additives. 



Except for blood plasma and some solutions used in 
variable pH studies, the fluids were bubbled with a 
5 % CO2 in air mixture to maintain the pH value at 
7,4. 

For dissolution-rate measurements, the cassettes 
were placed into a water bath at 37.0 **C ± 0.5 K, and 
fluid was pumped through them at a rate of 
(0.20 ± 0.04) ml/nun by peristaltic pumps. To avoid 
leakage, it was necessary to tap the cassettes for 
threaded fittings and to clamp the cassettes together. 
Except for short sections in the pumps, ail tubing was 
TFE (teflon) to avoid gas diffusion through the 
tubing. For selected samples, the fluid pH value was 
monitored by flow-through pH electrodes both 
before and after the cassettes. The solution pH value 
prior to the cassettes was 7.4 ± 0,2. Solution pH after 
it had passed through the cassettes was initially 
between 8 and 9, but for most samples it dropped to 
7.4 ± 0.2 within several hours. With the most rapidly 
dissolving samples it remained in the range 7.6 to 7.9 
for the rest of the experiment. This observation 
indicates that the flow rate used in this study may 
have been too low for the rapidly dissolving glasses 
and that for these glasses, dissolution rate may 
change if flow rates were incr ased. Further tests will 
be run to determine if this problem existed and its 
etfect, if any, on the measured fiber dissolution 
rate. 

Hie durations of th experiments were generally 
1, 2, 3, and 4 weeks. Individual samples were used for 
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Table Z Calculated 


glass fiber oxide compositions^) in wt% 
















glass 
no. 




AI2O3 


CaO 


MgO 


NazO 




B2O3 


Ti02 




SrO 


BaO 


F2 


1 


6132 


8.06 


9.18 


4.33 


15.31 


0.80 


0.50 


0.10 


0.30 


— 




— 


2 


54.00 


9.90 


15.80 


5.60 


9.20 


1.20 


4.30 


— 


— 


— 


— 


— 


3 


62.00 


8.00 


8.00 


4.00 


16.00 


— 


2.00 


— 




— 


— 


— 


4 


58.00 


8.00 


8.00 


4.00 


16.00 


— 


6.00 




— 




— 


*— 


5 


63.10 


5.50 


6.00 


3.10 


14.80 


1.10 


5.40 




0.10 


— 


— 


0.70 


6 


66.00 


4.00 


8.00 


4.00 


16.00 


— 


2.00 


— 


— 


— 


— 


— 


7 


64.15 


3.71 


7-98 


3,63 


15.28 


0.75 


4.43 


0.19 


0.30 


0.08 


— 


— 


S 


61.64 


4.00 


7.72 


3.72 


14.85 


0.84 


3.60 


0.27 


0,32 


0.09 


2.40 


— 


9 


59.51 


4.98 


736 


3.64 


14.22 


2.10 


7.53 


0.07 


0,37 


0.06 


— 


— 


10 


59,37 


5.14 


7.98 


3.82 


14.29 


2.17 


6.61 


0.08 


0,37 


0.04 




— 


11 


62.00 


4.00 


8.00 


2.00 


18.00 


— 


6.00 


— 


— 


— 


— 


— 


12 


62.00 


4.00 


8.00 


8.00 


.12.00 


— 


6.00 




— 


— 


— 


— 


13 


62.00 


4.00 


8.00 


4.00 


16.00 


— 


6.00 


— 


— 


— 


— 


— 


14 


59.20 


4.50 


15-30 


5.60 


10.90 


0.60 


3.90 


— 


— 


— 


— 


— 


15 


60172 


3.81 


7.92 


3.65 


14.78 


0.93 


5.07 


0.13 


0.26 


0.04 


2.02 


— 


16 


62.00 


4.00 


10.00 


2.00 


16.00 


— 


6.00 


— 


— 


— 


— 


— 


17 


60.50 


3.60 


15.90 


3.60 


10.50 


0.30 


5.60 


— 


— 


— 


— 


— 


18a 


61.99 


3.93 


7.91 


3.50 


14.63 


0.77 


6.97 


0.05 


0.17 


0.08 


— 


- 


18b 


61.99 


3.93 


7.91 


3.50 


14.63 


0.77 


6.97 


0.05 


0,17 


0.08 


— 


— 


18c 


61.99 


3.93 


7.91 


3^0 


14.63 


0.77 


6.97 


O.OS 


0.17 


0.08 


— 


— 


19 


62.00 


4.00 


4.00 


8.00 


16.00 




6.00 










— 


20 


60.16 


3.66 


8.09 


3,14 


14.42 


0.91 


6.04 


0.19 


0.25 


0.12 


2.45 




21 


62.00 


4.00 


12.00 


4.00 


12,00 




6.00 












22 


62.00 


4.00 


6.00 


4.00 


18.00 




6,00 












23 


64.71 


2.10 


8.60 


3.60 


14.90 


0.43 


5.40 


0.07 


0.12 


0.02 






24 


60.00 


4.00 


8.00 


4.00 


16.00 




8.00 












25 


59.39 


3.36 


7.97 


3.59 


14.23 


0.60 


9.78 


0.34 


0.32 








26 


64.00 


2.00 


8.00 


4.00 


16.00 




6.00 












27 


65.69 


0.50 


5.66 


3.02 


15.78 


0.80 


8.04 


0.10 


0.30 








28 


62.00 


2.00 


8.00 


4.00 


16.00 




8.00 












29 


66.00 




8.00 


4.00 


16.00 




6.00 












30 


62.94 


0.90 


8.70 


1.18 


16.06 


0.18 


9.56 


0.03 


0.18 


0,20 







^) Glass compositions are calculated from the raw material compositions with no attempt to correct for volatilization during melting. Glasses 
no. 18a to c are batches of filler made from the same bulk glass at different times. 



each time. At the end of the allotted time, the sample 
was removed from the water bath, cooled to room 
temperature and washed by switching the pump fluid 
to deionized water for at least 2 h. The cassette was 
then disconnected from the pumping lines, the water 
removed by suction filtration, and the cassette dried 
to constant weight in the manner already described. 
For several samples the dissolution was monitored for 
times up to 20 weeks to assure that the standard 4 
week experiments gave results representative of 
dissolution at longer times. 

SEM analyses of the partially-dissolved samples 
were perfonned on 1 cm diameter cores taken 
verticaUy at the center of each sample through the 
entire thickness. The core was cut in half lengthwise 
and embedded in epoxy resin. It was then ground 
briefly on the face of the cut and polished with 
diamond paste. This gives a polished section perpen- 
dicular to many fibers. 

Chemical analysis of dissolution fluid and glass 
fibers was done by inductively-coupled plasma arc 
spectroscopy (atomic absorption spectroscopy for 
potassium). Dissolution fluid samples were collected 
in 50 ml polyethylene bottles direcdy from the 



outflow tubing. They were immediately fix>zen and 
stored in the freezer until analyzed. Glass fibers were 
dissolved in a mixture of aqua regia and hydrofluoric 
acid in a microwave oven prior to analysis of all 
elements but boron. For boron analysis the boron was 
brought into solution by leaching the fibers in 
hydrofluoric acid. 

3. Dissolution rate constant calculation 

Glass fiber dissolution rates were calculated from 
weight-loss data assuming that the rate is propor- 
tional to the reaction-surface area: 

-dM/dt^kA (1) 

where M = fiber mass in ng at time r = time in h, 
k = dissolution-rate constant in ng/(cm- h), A = re- 
action-surface area in cm^. 

If the distribution of fiber diameters in the sample 
is assumed to be Gaussian, the following relation can 
be derived (section 8.1.): 

A£ ^ b exp (-^^) + (V?^/2) (H-2 b'^) (1+erf b) n) 
Afo ^0 exp i-bl) + (VS/2) (1+2 bl) (1+erf bp) 



18 



Gtastech. Ber. 64 (1991) Nr. 1 



Glass fiber dissolution in a physiological saPtne solution 



Table 3. Glass fiber ph^cal properties 

glass density^) fiber o in |un k in ng^ a in ng^ type of 
no. in g^cm^ diame- (cm^h) (cm^h) moTphoIogy 







ter 

in |un 








(see text) 


1 


2^2 


8.85 


0,84 


0.9 


0.2 


— 


2 


2^9 


9.40 


030 


3.0 


03 


— 


3 


2.50 


8.92 


031 


5.4 


1.4 


— 


4 


2.52 


939 


0.33 


9.0 


2.1 




5 


2.51 


9.62 


0.45 


27.4 


1-1 


A 


6 


2.49 


9.41 


0.39 


42.7 


3.4 




7 


2.51 


8.21 


032 


49.2 


4.4 


— 


8 


2^ 


9.48 


1.50 


57.5 


2.5 


— 


9 


2^1 


8.88 


1.02 


77.7 


6.8 


A 


10 


2.52 


9.15 


0.97 


84.1 


11.3 


^ 


11 


2.53 


9.25 


0,24 


108.0 


4.8 


A 


12 


2^2 


8.84 


0.42 


109.9 


11.9 


A 


13 


2.51 


9.97 


0.32 


139.0 


4.8 


A 


14 


2.59 


9.21 


0.39 


142.9 


12.3 


— 


15 


236 


9.27 


0.99 


150.1 


20.6 


— 


16 


233 


9.24 


0.28 


166.6 


113 


A 


17 


2.59 


9.56 


033 


202.2 


1.6 


A 


18a 


231 


8.26 


030 


204.7 


5.6 


A 


18b 


231 


9.08 


0.46 


85.6 


3.1 




18c 


231 


8.96 


0.23 


116.4 


2.1 


A 


19 


2.49 


9.24 


0.28 


206.0 


18.2 




20 


237 


9.73 


1.17 


214.1 


20.0 




21 


234 


9.27 


0.29 


218.8 


11.5 


A 


22 


2.51 


9.26 


038 


254.6 


26.1 


A 


23 


232 


9.67 


0.74 


2673 


26.4 


A+B 


24 


232 


9.33 


0.50 


375.9 


19.4 


A+B 


25 


233 


9.41 


0.75 


462.1 


34.0 




26 


2.51 


8.94 


038 


586.2 


36.8 


A+B 


27 


230 


8.80 


1.11 


596.5 


16,7 


A+B 


28 


233 


9.25 


031 


764.6 


68.1 


B 


29 


232 


9.27 


0.94 


765.9 


36.5 




30 


2.53 


9.08 


0.43 


886.9 


67.4 


B 



^) Fiber densities are calculated frozn the compositions in table 2 
based on a large number of measured values. Their zcaax^sy is 
±0.Q2g^cm^ 



where = initial fiber mass in ng, b = 
= {dp - 2 k //^)/(V2 o), bo = dol{\/2 o),do-- aver- 
age initial fiber diameter in jxm, g = fiber density in 
ng/cm^, a = standard deviation of the fiber diameter 
in pjn, erf J' = error function of y is 

erf y J exp (-jr?) dx. 

0 

The dissolution-rate constant, k, can be calculated by 
fitting the fiber weight loss data to equation (2) using 
non-linear least squares techniques. In section 4. it 
will be seen that 

a) the analysis gives a good fit to the data; 

b) the assumptions of the analysis are met for most 
glass fiber compositions despite a complex dissolution 
process; 

c) the analysis gives dissolution-rate constants that 
appear to describe the in vivo dissolution process. 

For the special case of a uniform fiber diameter, 
equation (2) reduces to: 

1 - iM/Mo)^ = 2 A: t/ido Q). (3) 




0 20 40 60 60 100 120 140 



1 


f 1 


1 


!■ 1 


0.04 


- b) 


fr=27 . 


qlass no.; ~ 


0.02 




^=0.9 


■ = 5 

a: 1 


0 — 


A 6 

1 1 







0 20 40 60 80 100 120 
f in d ^ 

Figures la and b. Dissolution data for various glass fibers. The lines 
are fit to the measured points by equation (2); the corresponding 
k values are given in ng^(cm^h), a) glasses no. 27, 18a and c, 
b) glasses no. 5 and 1. The point not used in the fit is enclosed in 
paraatheses. 

This is the equation describing dissolution of fibers at 
a constant dissolution velocity (constant thickness 
dissolved per time), which . follows directiy from 
equation (1) in the case of uniform fiber diameter 
[3 and 4]. Equation (3) is a close approximation to 
equation (2) for most of the sample s measure d here, 
so that plots of time versus 1 — V(Af /A^) are ap- 
proximately linear. Rate data are plotted in this way 
in this report but the dissolution-rate constants and 
the lines through the data have been calculated using 
equation (2). 

4. Dissolution mechanism 

Several fiber samples were examined in detail to 
assure that the method of data analysis is applicable 
to the dissolution process. The dissolution rates of the 
glasses no. 1, 5, 18a, c, and 27 span the range of 
dissolution rates measured (tables 2 and 3). TTie 
dissolution data for these glasses are fit well by 
equation (2) (figures la and b), but there is one effect 
which needs further comment. For each of the five 
glasses there is an initial period of slower dissolution 
rate. Scholze and Conradt [4] have noted a similar 
effect. Whatever the cause, it means that rate 
constants calculated from only the initial month of 
dissolution will generally underestimate the rate. For 
three of the glasses in figures la and b the initial rate 
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Figtues 2a to d. SEM backscattered electron images of glass no. 18a fiber after various dissolution times, /, showing unaltered glass, G, and 
leached layer, L. Morphology type A, a) untreated, b) I = 21 d, c) I = 55d, d) f = 84 d. 



is the same as that from longer times within the error 
of the measurement; for two, it is significantly smaller 
(glass no. 5: 21 %; glass no. 18c: 13 %). The error is 
not large compared to the variation with glass 
composition. The dissolution of most glass compo- 
sitions was therefore followed for only one month. 
This allowed more compositions to be examined so 
that a better understanding of the effect of glass 
composition could be developed. 

The partially-dissolved samples of the five glasses 
were examined under the SEM to correlate weight 
loss with morphological changes in the fibers. No 
changes were observed for glasses no. 1 and 5. 
Apparently the dissolution rates for these glasses are 
too slow to show morphologic effects within the 
approximately three-month period of measurement- 



Samples of glasses no. 18a and c (fiberized from the 
same glass sample) have the same changes in 
morphology: a leached layer separated from an 
unaltered core by a sharp boundary (morphology A, 
see figures 2a to d). The SEM backscattered electron 
images, in which brightness is dependent on the 
density and atomic mass of the material, indicate that 
the leached layers arc lower in density and/or average 
atomic mass than the core. Energy dispersive x-ray 
analysis confirms that concentrations of sodium, 
magnesium, and calcium are much reduced in the 
leached layer, leaving it enriched in silicon and 
aluminum. The SEM was unable to resolve finer 
variations in chemistry. The leached layers are all of 
approximately uniform thickness around the circum- 
ference of any given fiber (fiber no. 2) and there is 
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I 1 10 pm 

Figuze 3. SEM badcscattered electron image of glass no. 30 fiber after dissolution for 28 d showing unaltered glass, C, leached layer, and 
protective layer, P, associated iMth fibers in var^iig stages of dissolution. Morphology type B. 



little variation in thickness of the layers for different 
fibers in the same sample. This suggests that sample 
preparation and solution flow rate preclude stagnant 
areas in the sample. The leached layer thickness 
increases linearly with time (r^ = 0.9989) while the 
overall fiber diameter remains essentially unchanged. 
It appears that the dissolution rate of the leached 
layer is so slow as to be immeasurable with the 
techniques used here. However, one can see in 
figures 2a to d that the outer layer began to swell 
(most likely due to hydration) and separate from the 
fiber. Other SEM pictures show outer layers com- 
pletely separated; this suggests that the leached layer 
is mechanically weak and easily removed by physical 
processes. The rate of increase in leached layer 
thickness with time of glass no. 18a is 1.8 times 
greater than that of glass no. 18c. Their dissolu- 
tion-rate constants determined from weight loss are in 
the same proportion, indicating that for these two 
samples, the increase in leached layer thickness is 
proportional to the rate constant. Tliis large differ- 
ence in leaching rate between two samples with the 
same composition is discussed in section 4.2.6. From 
the data in figures la and b and 2a to d and the fiber 
geometry, the density of the leached layer can be 
calculated. It is approximately 0.5 g/cm^, corre- 
sponding to a loss of 80 % of its original mass. 

Some fibers of glass no. 27 have a dissolution 
morphology similar to that of glasses no. 18a and c. 



but others have a different dissolution morphology. 
This morphology (morphology B) is most clearly 
illustrated by fibers of a different composition (glass 
no. 30) which is shown in figure 3. What initially 
appears to be a bewildering variety of forms, can be 
understood in terms of a simple set of processes. The 
initial attack of the solution forms a thin leached layer 
of the order of 0.5 |xm in thickness. This layer protects 
the fiber from further dissolution so that in some 
instances the fiber remains unaltered, retaining its 
original diameter, even after months in solution. 
More commonly, the layer is breached in places and 
under the layer rapid dissolution occurs radially into 
the fiber ft-om these points. Thick leached layers are 
formed followed by complete dissolution of the fiber 
material under the outer protective layer. Eventually, 
the protective layer may disintegrate, leaving a small, 
leached remnant of the original fiber, which then 
dissolves completely. Energy dispersive x-ray analysts 
shows that the leached material is depleted in sodium, 
calcium, and magnesium. Aluminum is in too low a 
concentration to be determined. No difference in 
composition could be found between the protective 
leached layer and that inside the protective layer. 
Although the dissolution process seems fairly obvious 
from the SEM micrographs, the chemical basis for 
these processes are not yet understood. 

The dissolution solutions from those samples kept 
longest in the fluid were sampled every 2 to 7 d during 
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Table 4. 


Major component maj 


^ balance for 


glass no. 


18c fiber 


dissolution after 27.8 d in mgfg 


initial fiber 








measured 




calculated 








initial 


dissolved 


dissolved 


leached 


core 




fiber 






layer 






620 


64 


64^) 


49^) 


507 


AI2O3 


39 


0 




7^) 


32 


CaO 


79 


16 


14«) 


0^) 


65 


MgO 


35 


6 


6«) 


0^) 


29 


NazO 


146 




27*) 


0^) 


119 


K2O 


8 




1^) 


0^) 


7 


B2O3 


70 


14 


13^) 


0^) 


57 


total 


1000 


125*) 


125 


56 


819 ' 



*) Based on total fiber-weight loss assuming no hydration of 

leached layer after drying. 

^) Based on the measured amount dissolved. 

^) Assuming total icmoval from the measured leached layer and no 

removal from the core. 




Figure 4. Fiber-weight loss calculated from solution ana^ts (line) 
compared to measured weight loss (points) for glass no. 18c 
fiber. 



the dissolution for glasses no, 1, 18c, and 27. The 
solutions were subsequently analyzed for AI2O3, 
B2O3, C:aO, MgO, and SiOi, the major fiber 
components not originally present in high concen- 
tration in the dissolution fluid. It was necessary to 
correct the CaO and MgO analysis for the amounts 
originally in the solution. No AI2O3 was found in any 
solution sample. 

The solution-analysis data can be combined with 
weight loss and SEM analyses to give a consistent 
picture of the dissolution process for several of the 
glasses. The SEM energy dispersive x-ray analyses 
suggest that only Si02 and AI2O3 remain in the 
leached layers and that Si02 is partially leached. For 
glass no. 18c, the leached layer volume was calculated 
as a function of dissolution time from the regular 
dissolution geometry. From the mass of unaltered 
fiber corresponding to this volume, the concentration 
of CaO, MgO, and B2O3 in the leachate was 
computed assuming that these oxides were com- 
pletely removed from the leached layer. To compute 
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Figure 5. Variation of leached layer thickness at 28d with 
dissolution-rate constant calculated by weight loss and normalized 
by dividing by density. The dashed line shows the theoretical 
relation based on 80 % weight removal from leached layer, the 
solid line the calculated reduction in diameter for 100% re- 
moval. 



a mass balance it was assumed that Ka20 and K20 
were also completely leached from the leached layer. 
A comparison of these values for glass no. 18c after 
dissolution for 27.8 d to the measured values from 
solution analysis (table 4, the second and third 
columns, respectively) shows excellent agreement. 
The sum of the dissolved components is also in 
agreement with the mass lost from the fiber after 
27,8 d and over the entire time during which the 
solution was monitored (figure 4). 

It is not possible to construct tables similar to 
table 4 for the other two glasses for which the 
dissolution fluids were analyzed. The leached layer 
thickness was not measurable for glass no. 1, and 
glass no. 27 exhibits dissolution morphology B. But 
plots similar to figure 4 can be constructed for these 
glasses since the plots do not involve dissolution 
morphology. The agreement between weight loss and 
solution-analysis data for glass no. 1 is within the 
error of the measurements for such a slowly dis- 
solving glass. Weight loss calculated from solution 
analysis for glass no. 27 is approximately one half of 
the measured values, indicating some error in the 
data or the analyds. It is likely that the high solute 
concentrations resulting from the rapid dissolution of 
glass no. 27 caused precipitation in the dissolution 
fluids prior to analysis, but this has not yet been 
established. 

A survey of the dissolution morphologies of 
selected glasses in tables 2 and 3 was conducted to 
determine if the results of the detailed studies could 
be generalized to other IBber compositions. A smooth 
transition in morphology was found in going from 
slowly to rapidly dissolving glasses (figure 5). Glasses 
with slow dissolution rates dissolve with morphology 
A. As glasses are examined in order of increasing 
dissolution rate, samples begin to show a few fibers 
which ar dissolving with morphology B . The number 
of fibers dissolving with morphology B increases 
with dissolution rate until it becomes the dominant 
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morphology for glasses with rapid dissolution rates. 
In addition to pure morphologies A and/or B, fibers 
may also exhibit intermediate forms. The factors 
controlling the shift in morphology from sample to 
sample and within a given sample are not yet 
understood and are under investigation. 

When the leached layer thickness of a variety of 
fibers is plotted against dissolution-rate constant 
calculated from weight loss and normalized by 
dividing by density (figure 5), fibers with morphology 
A cluster around a line calculated assuming 80 % 
weight loss from the leached layer as was found for 
glasses no. 18a and c. Thus suggests that when 
morphology A is found, the dissolution process is 
basically the same regardless of the glass fiber 
composition. Glasses exhibiting morphology B, 
where there are some regions of complete dissolu- 
tion, he below this line and closer to the line 
describing the reduction in diameter brought about 
by complete dissolution. 

The mathematical model by which the dissolu- 
tion-rate constant is extracted from the weight-loss 
data is based on the assumption that the rate of 
weight loss is proportional to the reaction-surface 
area. The work on the dissolution mechanism 
supports the validity of this assumption for those glass 
compositions dissolving with morphology A. The 
theoretical fit to the weight-loss data is good, and the 
increase in leached layer thickness with time is 
predicted by . the model. This linear dependence 
eliminates the possibility of control by diffudon 
across the leached layer. The mathematical model is 
not completely correct, however, since it assumes 
complete rather than partial dissolution of the fiber. 
Equation (2) can be easily modified for the case of 
firstly, leaching to leave an insoluble layer or 
secondly, leaching to leave a layer which slowly 
dissolves (see section 8.)- There is uncertainty about 
whidi of these modifications, if any, to use for the 
following reasons: 

a) There remains some doubt about the fraction of 
material remaning in the leached layer and the effect 
of fiber composition and other variables on this 
fraction. 

b) There have been no measurements of the disso- 
lution rate of the leached layer. 

Figure 6 shows the fit of the extreme assumptions 
— equation (2) assuming 100% dissolution and 
equation (9) (section 8.2.) assimiing an insoluble 
fraction. The data and the curves have two basic 
regions — an early, approximately linear region and 
an approach to a constant level. The data collected 
for nearly all samples fall within the linear region 
where the two equations fit the data equally well. 
Glass no. 27, which is illustrated, has a large 
dissolution-rate constant and was followed for a long 
time. It exhibits both morphologies A and B. The last 
4 data points clearly require some insoluble or much 
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Figure 6. Comparative fits to weight4oss data for glass no. 27. The 
solid line assumes 100 % dissolution, the dashed line a 9 % 
insoluble fraction. 

more slowly dissolving fraction. The fit to these 
points is primarily a function of the firaction which is 
insoluble; the rate constant is the same for both lines. 
The scatter in these points for this sample may be 
related to varying proportions of the two morphol- 
o^es. The insoluble fi^ction fit to the data is an 
average of the last four points: 9 %. Glasses leaching 
with morphology A would be expected to have 
« 20 % residual material. Glass no. 18c was also 
followed for long time (figure la). It has a rate 
constant about V5 that of glass no. 27 and leaches with 
morphology A. Fitting the data to equation (9) 
(section 8.2.) with 20 % insoluble fraction gave a 
dissolution-rate constant 5 % greater than from 
equation (2). 

Most samples were not followed long enough to 
approach a constant level from which an insoluble 
fraction could be determined, nor is the fit to the 
availaSle data improved by using a more complicated 
equation. Because of these factors and the uncer- 
tainties already discussed, equation (2) was used to 
obtain all rate constants presented in this report. The 
neglect of the leached layer should introduce no more 
than a 5 % error. 

Equation (2) is not expected to describe the 
dissolution of glass fibers dissoMng with morphology 
B since the changes in surface area with time are not 
predictable. The fit to the data is good (figure 6), 
however, and in the absence of a better model, the 
dissolution-rate constants were calculated by means 
of equation (2) for these glasses also. 

4.1 . Comparison to in vivo data 

The literature contains many in vivo studies of glass 
fiber dissolution, but only three which used fibers of 
initially uniform diameter and measured fiber di- 
ameter change with time [11 to 13]. These provide the 
most accurate measurement of in vivo fiber dissolu- 
tion rate since they avoid uncertainties associated 
with a broad fiber diameter distribution. The com- 
position of the fibers used in all three studies is that of 

23 



Russell M. Potter; Stephanie M. Mattson: 



Tabic 5. Fiber dissolution rate constants for glass no. 5 calculated 
from literature in vivo data 



fiber animal fiber k in reference 

dimensions location ng^(cfn^ h) 

diameter x in lung 

length in \im 



1-5 X 5 


rat 


unspecified 


1.7 


[11] 


l^x 5 


rat 


tissue 


2.1 


[13] 


1^ x60 


rat 


unspecified 


8.7 


[111 


1^x60 


rat 


tissue 


4,1 


[13] 


1.5X60 


rat 


macrophages 


10,4 


(13] 


3.0 X 10 


rat 


unspecitied 


0.0 


[11] 


3.0 X 60 


rat 


unspecified 


17.4 


[111 


3.0 X 10 


hamster 


unspecified 


15.2 


(121 


3.0 X 30 


hamster 


unspecified 


48.5 


(121 


3.0x60 


hamster 


unspecified 


48.2 


[12] 



in vitro dissolution-rate constant, k, from this study: 
27.4ng^(cni2h) 



T 1 1 r 
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1/r in 10"*^ K*^ fc- 

Figure 7. Arrhentus plot for dissolution of glass no. 18b fiber 
showing the least-squares fit to the data. 



glass no. 5. Dissolution-rate constants were calcu- 
lated from the data in these studies assuming a linear 
decrease in fiber diameter with time. These calcu- 
lated dissolution-rate constants appear in table 5 
together with other pertinent experimental data from 
these studies and the in vitro dissolution-rate constant 
for glass no. S measured in this study. The in vivo 
dissolution rate varies greatly with local environment 
in the animal and with animal species. The in vitro 
rate is well within this range. This agreement is not 
proof that the in vitro rates are good measures of in 
vivo dissolution rates, but it is comforting that the 
agreement is as good as the quantity and quality of 
the data permit. Clearly, in vivo measurements of 
dilfferent compositions are necessary to evaluate fully 
the accuracy of the in vitro results. 

Fibers recovered from animals show a different 
morphology from those generated in in vitro ex- 
periments. There is a reduction in diameter with no 
evidence of a leached layer. Scholze [5] attributes tlus 
to the in vivo presence of organic complexing agents 
which remove the leached material as it is formed. 
Alternatively, its mechanical weakness may lead to 



physical removal while in the animal or during the 
extraction process. The absence of the leached layer 
suggests that the leaching process is the rate-deter- 
mining step in vivo as it is in this study in vitro. The 
dissolution-rate constants measured in vitro would 
then underestimate the in vivo rate in ng/(cm^ h) by 
about 30 % in the worst case with 20 % retained in 
the leached layer. 

4.2. Effects of experimental parameters 

Some work was devoted to the effects of several of 
the experimental parameters on the dissolution-rate 
constant in order to develop a rough feel for their 
importance. This section summarizes the results. 

4.2.1. Solution composition 

The physiological saline solution used routinely in 
this work is only an approximation to what is thought 
to exist in the lung. The dissolution rates of several 
samples were measured in two additional fluids to 
give an indication of the importance of fluid com- 
position. These fluids are human blood plasma and 
that developed and used by Scholze and Conradt '[4] 
with the addition of 0.1 % NaN3. The compositions of 
these fluids and the glass dissolution-rate constants 
measured in them are summarized in table 5. The 
data are limited but suggest that fluid composition 
can change the rate constant by at least a factor of 
two. In addition, some measurements were made in 
deionized water (2 glasses) and 0.1 mol [Tris/HCl]/l 
pH 7-4 buffer (9 glasses). The rate in deionized water 
is about V4 that in the modified Kanapilly solution; 
that in the Tris buffer, about Viq. The dissolution-rate 
constant for glass no. 18a measured in the Kanapilly 
solution without the NaN3 additive was within 
experimental error the same as that measured in the 
modified Kanapilly solution. 

4.2.2. Temperature 

Dissolution rates were measured routinely at 37 ^^C. 
Although higher temperatures offer the advantages 
of increased dissolution rate and reduced bacterial 
and algal growth, these advantages are outweighed by 
the uncertainty of extrapolation to body temperature. 
Also, the physiological saline solutions tend to be less 
stable toward precipitation at higher temperatures. 
The effect of temperature was measured in the range 
25 to eO'^C for glass no. 18b. The data appear linear 
in an Arrhenius plot (figure 7), yielding an activation 
energy of 64 kJ/moL 

4.2.3. Flow rate 

For accurate results from a flow-through technique, 
the solution around the fibers must be exchanged 
rapidly enough so that the concentration of materials 
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Figure 8. Effect of pH value on dissolution-rate constant, for 3 



Table 6. Comparison of dissolution-rate constants, k, for different 
fiber diameters of various glasses under different pH condi- 
tions 



glass no. 


pH value 


k in ng^{cm^ 


h) 


difference 
in % 






at I0|un 

fiber 

diameter 


at 2|im 

fiber 

diameter 


18a 


7.5 


204.7 


260.0 


27 


" \ 


' 6.3 
7.5 
. 8.1 


24.0 
139.0 
285.9 


30.7 
156.0 
308.6 


28 
12 
8 


- 1 


' 6.3 
7.5 
. 8.1 


274.8 
765.9 
950.1 


247.3 

639.0^ 

879.6') 


10 
17 
7 



Based on only one point. 



dissolved from the fibers does not influence the 
dissolution rate. Large flow rates with uniform 
bathing of fibers are thus desirable. Balancing this is 
the need to keep flow rates low to limit consumption 
of the dissolution fluids. The dissolution rate of glass 
no. 18a was measured for flow rates in the range 0.08 
to 1.00 ml/min. Within this range the dissolution rate 
is independent of flow rate indicating that these rates 
are sufficiently rapid to eliminate the effects of 
dissolved components for these glasses. The standard 
flow rate of 0.2 ml/min was chosen to be well above 
the lower limit of this range. The fact that the fliud 
was generally unchanged in pH value after passing 
through the sample cells confirm that this is a good 
choice of fiow rate for all but the most rapidly 
dissolving compositions. 

4.2.4. pH value 

The dissolution rates of three glass fiber compositions 
were studied as a function of the pH value in the 
range 5 to 8. The fiber compositions were chosen to 
include a rapidly dissolving glass, a slowly dissol^g 
glass, and one in between. The pH value of the fluid 
was changed with the minimum possible change in the 
fluid composition. Within limits it was possible to 
change the pH value of the standard fluid simply by 
changing the composition of the bubbling gas: 0.5 % 
CO2 in air for a pH value near 8 and 50 % CO2 in air 
for a pH value near 6. In order to decrease the pH 
value to near 5, NaHCOa (routinely 2268 mgfl) was 
eliminated from the fluid, and it was bubbled with 
pure CO2. Glass fiber dissolution rate was found to 
increase strongly with the pH value (figure 8 and 
table 6). 



4.2.5. niter-pore diameter 

Scholze and Conradt [4] reported that filter-pore 
diameter strongly affected their rates as measured by 
silica analysis. In the present experiments, reduction 



in vacuum-filtration rate after partial dissolution 
suggested that the filters were generally retaining 
partially-dissolved fine material, but SEM analysis 
never showed any build-up of fine material against 
the filter. Variation in filter-pore diameter from 0.1 
to 0.8 imi for glass no. 18b had no effect on the 
measured dissolution rate. 



4.2.6. Fiber diameter and fiber-forming cx)nditions 

Glass wool insulation is produced by several manu- 
facturing processes and has a broad distribution of 
ctiameters with a mean of ca. 5 (im. Ideally the fiber 
diameter for dissolution studies should be in the 
range trapped in the lung, less than 3 [un and have a 
narrow distribution. It is difficult and time-consuming 
to form fibers with a uniform diameter this small, and 
they are difficult to disperse adequately in the sample 
•cells. Therefore, fibers approximately lOjtm in 
diameter were used routinely. If the assumptions 
underlying equations (1 and 2) are correct, dissolu- 
tion-rate constants calculated from equation (2) 
should be the same regardless of the diameter of the 
fiber used in the experiment unless differences in 
forming conditions change fiber properties affecting 
dissolution rate. For instance, the more rapid quench- 
ing rate associated with smaller diameter fibers might 
lead to a more open glass structure, causing an 
increase in dissolution rate. To evaluate such effects, 
several glass compositions were formed into fibers 
near 2|Am in diameter and the dissolution-rate 
constants measured. The rates are compared to those 
from the 10 |jim diameter fibers in table 6 for several 
different pH values. Glasses no. 13 and 18a have 
consistently higher rates for the 2fmi fibers. The 
reverse is true for glass no. 29. These differences are 
within the 95 % confidence limits of the dissolu- 
tion-rate constants, but the systematic differences 
seen for each glass type suggest that the differences 
are real; Hie dissolution rates of the 2 \im fibers of 
glass no. 29 are probably anomalously low due to 
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poor dispersion of the fiber in the sample. Poor 
dispersion has a greater effect on the more rapidly 
dissolving glasses, and SBM analysis suggests that this 
has occurred for the 2 |xm fibers of glass no. 29. 
Further work should be done on the effect of fiber 
diameter^ but for the present purposes, the agree- 
ment in dissolution rate for different diameter 
samples shown in table 6 is sufficient to allow use of 
10 ym fibers for routine work. 

The various processes used to manufacture glass 
wool insulation enforce different forming conditions 
such as quench rate. Although these results on 
fiber-diameter effects suggest that the fiberKiissolu- 
tion rate is not strongly influenced by forming 
conditions, there is some evidence to the contrary. 
Throughout this study, glasses no. 18a to c were used 
as standards to measure the effects of various 
xperimentai parameters. Twice the supply of fibers 
was exhausted, and it was necessary to form more. 
Thus, there is data on three fiber batches of glass 
no. 18, all made from the same bulk glass sample. 
The dissolution rates measured on these samples vary 
by over a foctor of 2 (table 3). This variation reflects 
differences in the fiber samples rather than in the 
experimental conditions since several other fiber 
samples were remeasured occasionally (glasses 
no. 13, 18a, 29) and their dissolution rates did not 
change. Chemical analysis of the three glass no. 18 
fiber samples showed no difference in composition, 
so the variation in dissolution rate must reflect 
differences in the physical state of the fiber. Work is 
in progress to understand this variation in rate. 



5. Effect of glass connposition 

The results described thus far suggest that the 
dissolution-rate constants calculated from the dissor 
lution experiments describe the in vitro dissolution of 
glass fibers with reasonable accuracy up to the point 
at which the entire fiber has been leached. A 
significant, unexplained variation in dissolution rate 
has been found that appears to result from differences 
in the physical state of the fiber, but this variation is 
still much less than the variation in rate across the 
compositional range examined (table 3). The tech- 
nique thus appears applicable for survey work on a 
large number of glasses to determine how the in vitro 
dissolution rate varies with composition. 

Tables 2 and 3 list the chemical compositions and 
dissolution-rate constants of the glasses used in this 
study in order of increasing dissolution rate. It is clear 
that chemical composition can have a profound effect 
on dissolution rate. Glass fibers which have nearly 
identical shape and size and which have been formed 
in the same manner, show a variation in dissolution 
rate of nearly three orders of magnitude. This 
variation must be ascribed to the compositional 
variation of the glasses, which is within a fairly 



narrow range in terms of the possible compositional 
variation of glass fibers. 

Some of the effects of the individual components 
are obvious from a visual inspection of the tables. The 
primary compositional effect appears to be that of 
AI2O3, which strongly reduces dissolution rate, as 
seen by comparing glasses no. 4, 13, and 26. This 
conclusion supports that of Bauer et al. [6]. In 
addition, B2O3 appears to increase dissolution rate 
(compare glasses no. 6 and 13), although the influ- 
ence is not so pronounced as that of AI2O3. A more 
systematic analysis by standard linear least-squares 
fitting techniques showed the influences of the other 
oxides. Although some complications arise due to the 
fact that the dissolution mechanism is gradually 
changing as the dissolution rate increases, the major 
oxide components clearly fall into three catego- 
ries: 

first, AI2O3: which strongly decreases dissolution rate 
by a factor approximately two times greater than 
those of the second group; 

second, B2O3, BaO, Na20, CaO, and MgO: which all 
increase dissolution rate to about the same degree, 
although the effect of B2O3 is slightly larger; 
third, Si02: which does not strongly influence the 
dissolution rate. 

Typical glass wool insulation is based on small 
variations around the composition given by Mohr and 
Rowe [14]. This compositional range is much nar- 
rower than that covered in this study. The dissolution 
rate of glasses within this narrow range was computed 
to vary between 50 and 300 ng^(cm^ h). 

As a comparison, the durability of chrysotile 
asbestos from the Thetford Mine in Quebec (Canada) 
was measured. This weight-loss data was analyzed 
using a modification of equation (2) for a single fiber 
diameter (0.084 |jun) estimated from surface-area 
measurements. The measured rate constant of 
(0.09 ± 0,05) is in good agreement with 
0.05 ng/(cm^ h) calculated ft-om silica loss in a 
different modification of Gamble's solution pres- 
ented by Scholze and Conradt [4]. It is much larger 
than 0.004 ng/(cm^h) — a value calculated from 
surface measurements by Bauer et al. [6]. All these 
values are much lower than that measured for glass 
wool. 

To see these glasses and the asbestos from 
another perspective, one can compute how long it 
would take a 1 M-m fiber to dissolve in simulated lung 
fluid at 37 ®C. For glasses with dissolution constants 
of 50 and 300 ng/(cm^ h), the time to total solution 
would be 83 and 15 d, respectively. In contrast, 
chrysotile with a dissolution-rate constant of 
0.09 ng/(cm^ h) would require 47 080 d or 129 years to 
dissolve. These computations assume that the leach- 
ed silica-alumina remnant observed in some in vitro 
studies is eliminated in a secondary process. This 
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assumption is supported by the absence of such a 
leached layer on fibers collected from lungs of 
animals or humans. 

6. Concluding remarks 

This work has examined the influence of chemical 
composition on the dissolution of glass fibers under 
conditions approximating those thought to exist in the 
lung. The dissolution rate of compositions dose to 
commercial glass wool varied by 3 orders of magni- 
tude. Several key effects are apparent: 

a) The strongest compositional effect is that of AI2O3, 
which gready retards dissolution. This is a 
well-known effect for dissolution in water and in 
neutral salt solutions [15] and has been indicated by 
earlier studies of glass fiber in physiological saline 
solution [6]. Ai203 is essentially insoluble in water 
solutions at a near neutral pH value, but it has a 
greater effect than simply failing to dissolve. It 
somehow retards the leaching process for other fiber 
components as welL 

b) Dissolved alkali-metal cations (excluding lithium) 
present in near-neutral solution accelerate glass 
dissolution. This effect has been known for some time 
[15] and has previously been noted for glass fiber 
dissolution in saline solution relative to distilled water 
[3]. The work presented here shows the effect to be 
even greater when the saline solution is compared to 
a Tris/HQ buffer, in which the pH value is main- 
tained at 7-4. 

c) Dissolution morphology has been observed to 
change as a function of dissolution rate. Neither the 
cause nor the mechanism is yet known, but the 
change in mechanism involves the formation of a 
protective layer which retards the dissolution of the 
more rapidly dissolving glasses. 

The applicability of these results to glass fiber 
dissolution in the lung remains uncertain. For the 
limited data available currentiy, there is good agree- 
ment between the measured dissolution rates in vitro 
and in vivo, and differences in dissolution mor- 
phology can be accounted for plausibly. There is a 
pressing need for quantitative measurements of in 
vivo fiber dissolution of a number of different fiber 
compositions so that it can be determined if large 
differences in in vitro dissolution rate do indeed 
correlate to similar differences in in vivo dissolution 
rate. Such comparisons may support the interpreta- 
tions of chemical and physical processes that are 
crucial to understanding dissolution of glass in 
vivo. 

The authors gratefully adcnowledge the following colleagues, who 
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Gilbert, P. A. Tavcner, T. L. Gano, and H. Spencer. 
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8. Appendix: Mathematical derivations for the calculation of 
fiber dissolution and weight loss 

8.1. General derivation 

This appendix contains the derivation of the relation between 
sample mass and time used to calculate fiber dissolution rate 
constants in this report (equation (2)). It also contains modifica- 
tions of this equation for two different dissolution conditions. 

If the loss of material from a given fiber is proportional to its 
surface area, then: 

-W = *^ " -|(-fa^e/.)=*»rfZ.. (4) 

In this and the following equations, k = dissolution-rate constant, 
M — fiber mass, A = fiber-surface area, d « fiber diameter, g = 
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fiber density, L = fiber length and / = time. When the term in 
parentheses is differentiated, neglecting the change of length with 
time, equation (4) reduces to: 



8/ g 



(5) 



which can be integrated to yield: 

d = do- 2k tig (6) 

where 4o = initial fiber diameter. This equation can be applied to 
ail fibers. 

The initial diameter of the fibers is assumed to follow a 
Gaussian distribution of mean and standard deviation a. If the 
total number of fibers is W, the probability of having a fiber of 
diameter do is 

A*) = cxp »/2 Udo - do)/a)^) . 

A mormalizing factor a has been introduced to take into account 
the fact that only positive values of d^ are being considered: 

1 * - 
« = 7^ / exp (- Vi {(rfo - do)/ay) Bdo . 

0 

The mass of a single fiber is jt^^ L/4, Substituting for d from 
equation (6) gives: 



After integration and simplification, equation (7) becomes: 

M _ b cxp (- b^) + (1 -H 2 b^) (1 + erf 6) 

Mo " 6oexp (- bi) + (\^/2) (1 +2bl) (1 + erf 6o) 

which is equation (2) in section 3. with 

6 = (^0 - 2kt/Q)/(y/2o); bo = dof(y/2a) and 



(8) 



2 f 

crf3^-^J cxp(-j:2)^ 



8.2. Calculation of fiber weight 
dissoluticxi 



loss curve for an iiKxvnplete 



In practice, for type A morphology* the fibers do not dissolve 
complcidy but a very porous structure containing alumina and 
silica remains. This structure dissolves very slowly. After complete 
leaching, a fraction fi of the initial mass Mq remains. The 
computation of the change in mass of the fiber with time can be 
simply treated by considering that the mass of the fiber is the sum 
of an insoluble part, ^ Mq, that remains constant with time and a 
soluble part, (1 - P) Mq, that dissolves according to equation (8). 
In this case, equation (8) becomes: 

M~0Mo b'cxp^ 
(1— /S)Afo ^ocxp 
where 



> (- b'^) A- (yj^fl) (1 -hlb'^) (1 -t- erf 6') 
p (- bl) + (V5/2) (1 + 2 6J) (1+ erf bo) 



The total mass of the sample at time t is obtained by integrating 
over ail diameters for </>0. From equation (6), this occurs for 
do>2k tig. Assuming that all fibers have the same length L: 

IktlQ 

Using the notation 

e(A) = 6e-** + (Vjt/2) (1 + 2 6) (1-1- erf 6). 
For the fast dissolving fraction, the mass Mx varies as: 
A#,/((I -fi)h4^^ Qibi)/Q(bo) 

where 6j = [Jo - 2 t/ig (I - /3)]/(V2<7) , (10) 
For the slow dissolving fraction, the mass M2 varies as: 
Af2/(|8i»^o) = G(*2)/fi{*o) 

where 62 = [rfo - 2 t/ig (1 - «l/(\^ cr) . (11) 
The overall equation is 

(12) 



6' = [rfo - 2 A: t/ig (1 - fi)]/iy/2 a) . 

In the following case, the leached layer described in section 3. is 
assumed to dissolve slowly. The fiber can be considered to consist 
of a rapidly dissolving fraction of iniUal mass (1 — fi) Afo, density 
(1 — fi) g, dissoWing at rate A:, and a slowly dissolving fraction of 
initial mass 0 Mq, density fig^ dissolving at rate /tj. For each 
fraction, equation (8) applies. 



(7) 



where the total mass W = Af 1 + A/2- *i <»n be computed at short 
times. k2 at long times. 



8.3. Derivation for the case of narrow fiber-diameter distribu- 
tion 

If 0, 6 « (unless t = g do/ (2 k)) and equation (8) simplifies 
to 



MlMo-^ibibof. 



(13) 



Substituting the definitions of b and_^o from equation (8) into 
equation (13), and recognizing that = do'. 



{M/Mfd^ ^ 1 - 2 * tf(do g) 
which is equation (3) in section 3. 



(14) 
91RO0Q2 
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